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Abstract. Recent developments in synthetic and supramolecular chemistry have 
created opportunities to design organic systems with tailored nanoscale structure for 
various technological applications. A key application area is the capture of light 
energy and its conversion into electrochemical or chemical forms for photovoltaic or 
sensing applications. In this work we consider cylindrical assemblies of chromophores 
that model structures produced by several supramolecular techniques. Our study is 
especially guided by the versatile structures produced by virus-templated assembly. 
We use a mult i- objective optimization framework to determine design principles and 
limitations in light harvesting performance for such assemblies, both in the presence 
and absence of disorder. We identify a fundamental trade-off in cylindrical assemblies 
that is encountered when attempting to maximize both efficiency of energy transfer 
and absorption bandwidth. We also rationalize the optimal design strategies and 
provide explanations for why various structures provide optimal performance. Most 
importantly, we find that the optimal design strategies depend on the amount of 
energetic and structural disorder in the system. The aim of these studies is to develop 
a program of quantum-informed rational design for construction of organic assemblies 
that have the same degree of tailored nanoscale structure as biological photosynthetic 
light harvesting complexes, and also have the potential to reproduce their remarkable 
light harvesting performance. 
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Interest in the molecular mechanisms underlying photosynthetic light harvesting 
has recently escalated, fueled in part by the potential of engineering hiomimetic solar 
energy harvesting technologies. The hiomimetic approach aims to reproduce properties 
of the light harvesting complexes (LHCs) found in biology by using solid-state or 
organic components engineered at the nanoscale. Natural LHCs are remarkably efficient 
at all of the primary stages of photosynthesis: light capture, energy transfer, free 
carrier generation, and charge transfer [Ij. Furthermore, LHCs perform these tasks 
in a manner that is robust to varying external and internal conditions. Reproducing 
such efficiencies and robustness would revolutionize our energy production capabilities. 
This constitutes part of the tremendous appeal of biomimetic approaches to designing 
photovoltaic technologies. However, photosynthetic light harvesting also has useful 
lessons for development of other technologies. In particular, efficient photon capture 
and conversion to charge is also an important component of sensor technologies, raising 
the question whether these might also benefit from biomimetic approaches. 

Biomimetic light harvesting may be viewed as one component, the "front end", 
of artificial photosynthesis, which seeks to generate energy rich materials or fuels from 
sunlight. Recent years have seen impressive achievements in several aspects of artificial 
photosynthesis, with progress in a diverse range of platforms ranging from molecular to 
semiconductor systems [2j. In this work we focus on light harvesting, the initial stage of 
any photosynthetic unit, and address the question of how the collection of light energy 
may be optimized by a biomimetic LHC and what tradeoffs are involved in achieving 
this goal, within a program of quantum-informed rational design. 

In order to construct a functional LHC from the bottom-up we require a detailed 
mapping between structural features and motifs, and the light harvesting or sensing 
capabilities of the composite system. Establishing such a mapping is complicated by 
the wide variety of different structures evidenced by natural LHCS. Despite the detailed 
variations in structure, however, some general motifs do emerge. In particular, the 
majority of LHCs consist of densely packed pigment-protein complexes in which the 
light-absorbing chromophores {e.g. chlorophyll and carotenoid molecules) are arranged 
with specific relative orientations and locations within scaffolds provided by proteins 
[H This complexity indicates that a sophisticated understanding of how the nanoscale 
structure infiuences the light harvesting function is necessary in order to introduce 
precision and accuracy into the biomimetic approach. In particular, recent studies 
of the behavior of the natural LHCs consisting of pigment-protein complexes have 
revealed evidence of quantum dynamical effects in the electronic energy transfer through 
the complex, indicating that quantum effects hitherto neglected in analysis of energy 
transport may play a role in the high quantum efficiency and need to be accounted for 

X Exceptions to this general feature are provided by the antenna complexes of photosynthetic bacteria 
that are adapted to survive under conditions of very weak illumination, e.g., green sulfur bacteria, 
which consist of very large numbers (10^ — 10^) of chromophores without a protein scaffold [3||4^. 
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[5]. Theoretical modeling of the optical and electronic properties of large-scale pigment 
protein structures is a key component of obtaining this mapping and building a deeper 
understanding of the structure-function relationships in light harvesting. 

Much recent theoretical analysis has focused on understanding the quantum 
efficiency of energy transfer in biological LHCs - e.g. Refs. [iSJ El El [HI [IDl [HI [12]. 
Detailed modeling, including incorporation of quantum coherent effects such as exciton 
delocalization and chromophore-protein interactions, has allowed rationalization of 
the energy transfer times seen in experiments and the high quantum efficiencies 
typical of LHCs. Theoretical simulations have identified the dynamical parameters 
necessary for optimal energy transfer in model LHC systems and shown that a delicate 
balance between quantum coherent and incoherent dynamics appears to characterize 
this optimal energy transport regime [131 [Ml- Calculations exploring the landscape 
of these parameters for a small, well characterized LHC, the Fenna-Matthews-Olson 
(FMO) complex, have further indicated that this component of the green sulfur 
bacteria photosynthetic system operates in just such an optimal regime [TOl [12] • This 
would suggest that the design of optimized biomimetic LHCs should engineer the 
nanoscale structure so as to achieve strong inter-pigment couplings that compete with 
decoherent processes such as vibrational relaxation. Since strong interpigment coupling 
is an important requirement for quantum coherence in electronic energy transfer, a 
related underlying question for the biomimetic LHC program is whether such quantum 
coherence is coincidental or whether it is essential to the light harvesting function. 

In this paper we take a different approach to the study of optimality in light 
harvesting by incorporating an important perspective from the field of multi-objective 
optimization. Our starting point is the recognition that light harvesting, whether by 
natural or artificial systems, is not uniquely focused on achieving a high quantum 
efficiency for conversion of light to charge carriers, but is inherently a multi-objective 
optimization with several key objectives that must be simultaneously taken into account. 
In biological systems there are a wide array of objectives, not all of which necessarily 
have the same weight or status. One could say that the most important objective 
for a biological system is survival, which will ultimately dictate the changes in design 
features in response to environmental changes. In constructing biomimetic LHCs, we 
have a simpler task in that while the number of objectives are still greater than one, 
they are probably fewer and more equivalent in rank than in the biological case. For 
example, for photovoltaic applications one will likely want not only to maximize energy 
transfer and free-carrier generation efficiencies, but also to maximize the spectral width 
of absorption. In contrast, for a sensor the second objective may be to maximize the 
sensitivity to a particular wavelength rather than the spectral width of absorption. It 
is well known that in the presence of such multiple objectives one can have competition 
between them, which results in the emergence of families of optimal solutions that 
negotiate the trade-offs between the competing objectives in different ways [15j. 
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The questions posed in this work are thus two- fold. First, are there fundamental 
trade-offs involved in light harvesting? Second, if so, how can these be negotiated by 
engineering the structure of light harvesting complexes? We address these questions here 
by explicit calculation of the simultaneous optimization of the two desiderata mentioned 
above for a light harvesting system that might be used for energy conversion; we seek to 
simultaneously optimize the efficiency of excitonic transport and the spectral width of 
absorption. We do this by studying a prototypical biomimetic light harvesting antenna 
and requiring that it absorb photons in as wide a spectral window as possible while also 
efficiently transporting the resulting excitation energy to regions of charge separation. 
The ultimate limits of spectral width are dictated by choice of pigments; that is, pigment 
transition energies largely determine the width of absorption profiles. Therefore, one 
method for increasing the absorption profile width is to include pigments with as 
many different transition energies as possible. However, this creates an energetically 
disordered aggregate with typically reduced transport efficiency. Thus there is a trade- 
offs or competition, between the two objectives. One way to negotiate this conflict is 
to use molecular aggregation; by employing strong Coulombic coupling, aggregates of 
pigments can broaden or sharpen absorption profiles as exemplified by the classic H- and 
J-aggregates absorption profiles [16j. In this work we shall explore the extent to which 
this provides an effective technique for negotiating the trade-off between efficiency and 
spectral-width. 

We pause to mention some previous work examining the optimization and design 
of molecular structure in the context of light harvesting. Fetisova conducted some early 
and far-sighted studies into the relationship between structure and function in light 
harvesting complexes and strategies for optimization of structure, e.g. [171 [iHl [I9] . Much 
of this work was conducted well before experimental evidence for dynamical quantum 
coherence and therefore considered classical models of light harvesting dynamics. More 
recently, Fingerhut et al have employed structural optimization to design synthetic 
centers for ultra-fast charge separation for artificial photosynthesis [20j. The charge 
separation dynamics was modeled classically using Marcus theory, and notably, they 
examined a multi-objective optimization landscape and considered trade-offs between 
the quantum efficiency of charge separation and other objectives. Knoester and co- 
workers have examined the excitonic and optical properties of general cylindrical 
aggregates of identical chromophores (with and without disorder) in several works, 
including Refs. [211 [22]. Finally Noy et al [23], and more recently Scholes et al [24j, 
have compiled summaries of insights gained from studying natural LHCs and discussed 
how they aid the design of artificial light harvesting systems. 

The remainder of this paper is structured as follows. Section [T] introduces the 
physical systems forming the cylindrical assemblies of chromophores that constitute 
the focus of our study. Section [2] then outlines our theoretical models of the structure 
and electronic excitation dynamics in these systems. Section |3] presents the results of 
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our multi-objective optimization studies in terms of achievable objectives and tradeoffs 
encountered. This is followed in section |4] by a detailed analysis of the structural and 
excitonic properties of the chromophore assemblies that optimize the objectives. Then 
in section [5] the preceding study and analysis is condensed into a set of design principles 
for light harvesting complexes, particularly cylindrical antennas. Finally, we conclude 
with a discussion and an outline of future work in section |6l 

1. Cylindrical chromophore assemblies 

We focus in this work on cylindrical molecular assemblies as a prototypical architecture 
for biomimetic LHCs. Many systems self- assemble or can be templated into cylindrical 
structures and several of these have been studied as candidates for artificial light 
harvesting systems. Examples include carbocyanine molecules with hydrophobic and 
hydrophilic side groups and porphyrin derivatives that aggregate into cylindrical 
structures fSSlES]. There are also examples of cylindrical molecular aggregates in natural 
LHCs, the most prominent one being the chlorosome complexes of green sulfur bacteria 

One particularly promising realization of a cylindrical aggregate is not formed 
from the direct aggregation of chromophores, but rather, from virus-templated assembly. 
Virus-templated chromophore assemblies are supramolecular complexes constructed by 
attaching chromophores to protein coats of viruses that then self- assemble into large, 
regular structures. The self-assembling protein coats are used as rigid scaffolds that 
guide the synthetic organization of chromophores. Structures templated using the 
tobacco-mosaic, M13 and MS2 viruses, amongst others, have been demonstrated for 
use in light harvesting [271 EHl ES] , drug delivery [30] , and battery technologies [3ll [32] . 
This method of templated self-assembly presents a particularly promising path towards 
producing engineered assemblies of chromophores with controlled nanoscale structure. 
The templated assembly allows for a greater degree of customization than direct 
aggregation of chromophores because the nanoscale structure can be controlled by the 
choice of chromophores, templating protein, and of chromophore-protein linker groups. 
Templated assembly has the further advantage that pigments which do not naturally 
self-aggregate can be used. Also, it should be noted that the virus-templated assembly 
process creates a protein-pigment structure and not a direct molecular aggregate; as 
noted above, the former are more common in natural LHCs. 

We consider here templated assemblies based on use of the tobacco mosaic virus 
(TMV) as a scaflFold. The protein of this virus can self-assemble into stacked disks 
or cylinders, depending on assembly conditions such as solution pH [33j. From a 
practical stand-point TMV assemblies are advantages because the assembly process is 
very well known and large quantities can be reliably produced. Various laboratories have 
demonstrated covalent attachment of pigments to various sites on mutated TMV protein 
monomers [271 [28] (see Fig. [l]), and hence shown that the inter-pigment distances can 
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be tuned by choice of attachment site. This is advantageous for studying the impact 
of quantum coherent effects on hght harvesting, because coherent dynamics wiU be 
more prevalent in densely packed structures where inter-pigment couplings are strong. 
Thus the ability to vary the inter-pigment distance controUably over a finely tuned and 
large range of values, with other features of the protein scaffold held constant, allows 
the synthesis of chromophore arrays with and without the potential for such coherent 
dynamics. 



a 



TMV monomer with 
mutation at site 123 



TMV monomer with 
attached chromophore 
at site 123 




Assembled TMV disk 
that templates circular 
array of chromophores 



(a) Templated assembly using TMV 




(b) Variety of sites on a TMV protein monomer 
at which mutations can be introduced to 
facilitate covalent attachment of chromophores 

Figure 1. The self- assembling TMV protein provides a scaffold for producing an array 
of chromophores with well defined inter-chromophore distances. 



2. Structural and dynamical model 

2.1. Structural model 

Our multi-objective optimization studies are made for an idealized but versatile model 
for cylindrical arrays of chromophores that is shown in Fig. [2j The model consists of 
disks, each with M chromophores attached at specific sites, as determined by an 
implicit protein scaffold or template. Chromophore location on sequential disks may 
be ofl[-set by a variable amount and the chromophores have variable orientation but 
they are restricted to all have the same orientation relative to their disk. This model 
allows us to simulate TMV-templated chromophore assemblies in cylindrical and disk- 
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like structures, which may have vertical or helical stacking of chromophores along the 
cylinder. The cylindrical array of chromophores will function as an antenna absorbing 
photons and transporting the resulting photo-excitation to one of its ends. This end of 
the cylinder (referred to as the "bottom" of the cylindrical aggregate) interfaces with 
a surface or electron acceptor and charge separation occurs at this interface. For the 
TMV-templated structures a distinct chromophore whose energy levels match well with 
the remainder of the chromophores and the surface work function can be attached at the 
end to facilitate the electron transfer event. In this first work, we will not be concerned 
with optimizing these subsequent charge separation and transfer events but rather with 
the optimization of the preceding energy transfer to the separating interface. 

In the present calculations we assume that the following degrees of freedom (design 
parameters) can be tuned: (i) the transition energy of the chromophores on each disk 
(i.e., each of the disks has a distinct energy, but the M chromophores on any single 
disk have the same transition energy), (ii) all chromophores attach to the TMV protein 
scaflFold with the same orientation and this orientation defines a transition dipole that 
is specified by two angles, 9 a tangential angle and cp a radial angle, (iii) there can be 
a degree of misalignment between neighboring disks and this defines a helical angle, 9h^ 
which allows introduction of a helical twist along the cylinder coordinate. See Fig. [2] for 
representations of these angles. This results in a total of + 3 controllable degrees of 
freedom, which may be used to produce a wide range of chromophore array structures 
that range from complete alignment along any axis, to helical arrays with variable 
numbers of helical strands, e.g., M = 1 with 9h > corresponds to a single helical 
strand, M = 2 with > to a double helix, etc. TMV can template a 1-helix at the 
appropriate conditions, however we will not study this structure in this work, and will 
instead focus on M = 17 as specified below. In addition, we note that the chlorosome of 
green sulfur bacteria is assumed to have a structure that is an M-helix [4] and therefore 
this model is fiexible enough to capture artificial and natural light harvesting structures. 

The calculations presented here make explicit comparison of two TMV-templated 
aggregates based on the stacked disk morphology, which are designed to have very 
different inter-chromophore distances. Upon aggregation each TMV "disk" consists of 17 
protein monomers. Assuming complete functionalization (attachment of chromophores) 
this results in M = 17 chromophores per disk. Each chromophore is attached to the disk 
at a specific location, which is determined by its binding site on a protein monomer: this 
can be varied by prior treatment of the protein, allowing variable attachment location 
and hence variable inter-chromophore distance, see Fig. [T] We compare here TMV103 
and TMV123, where the numerical labels refer to the protein monomer site where 
pigments are attached. The 103 site is within the pore of the assembled TMV protein 
and attachment here results in a disk of chromophores of radius 25 A, implying an 
average neighboring pigment separation distance (within a disk) of about 10 A, and 
thus strong coupling of chromophores. In contrast, attachment at site 123 results in 
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Figure 2. Reduced description of chromophore assembly resulting from stacked disk 
TMV structure, in which only the chromophore properties are variable and the TMV 
protein scaffold is held constant. There are M chromophores attached to each disk and 
N disks. All chromophores are assumed to be attached with the same orientation and 
this orientation defines a transition dipole for each chromophore which is specified by 
two angles, 6 and (j). Finally, an angle Oh specifies the degree of misalignment between 
neighboring disks, and allows for a helical twist. 



a disk of chromophores of radius 40 A and an average neighboring pigment separation 
distance of 14 — 15 A, which leads to weaker electronic couplings. The vertical distance 
between disks remains the same for both structures: ~ 20 A. While TMV103 can achieve 
a greater density of pigments, we note, however, that the coupling between pigments 
on neighboring disks can nevertheless still be the same for both structures because 
the distance between disks is the same. A question we seek to answer in the following 
analysis is whether the dense packing within disks and any resultant increase in coherent 
dynamics of TMV103 is advantageous for any aspect of light harvesting. 

We note that current experimental techniques do not allow independent tuning of all 
the above design parameters for TMV-templated cylindrical assemblies. For instance, it 
is difficult to restrict a disk to only containing chromophores of a given species (transition 
energy) when there are several species on the whole cylinder. However, it is important 
to choose a large number of independent degrees of freedom in order to explore the 
landscape of optimal light harvesting given near-ideal resources. The figures of merit 
obtained can therefore be considered upper bounds on what is currently experimentally 
achievable. In a forthcoming publication we will examine the optimization landscape 
for a more modest set of design parameters. 
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2.2. Dynamical model 

The TMV-templated cylindrical chromophore assemblies are complex biomolecular 
systems. In order to efficiently compute and compare their light harvesting properties 
we must resort to an effective description of the electronic degrees of freedom. In analogy 
with the conventional modeling of natural light harvesting complexes jH we employ 
the Frenkel Hamiltonian to describe coherent dynamics: 

H^ = Y.mi^\ + Y,j,,{\^{j\ + \j){t\) (1) 

where \i) denotes an electronic excitation localized on pigment i in the complex, Ei is 
the transition energy of pigment i and Jij denotes the Coulombic coupling of pigments i 
and j. For simplicity, we assume that each pigment has only one dominant transition in 
the wavelength region of interest. This model can be easily generalized to pigments with 
more than one transition. In this work we treat the electronic coupling in the dipole- 
dipole coupling approximation. This approximation is generally valid for inter-pigment 
separations greater than ^ 12 A [35l [36], and while we will model some cases where the 
inter-pigment distance is slightly smaller than this, the error due to the dipole-dipole 
approximation in those cases will not affect our conclusions significantly. The eigenstates 
of i/g are the excitonic states in the complex, whose decomposition in terms of the site 
basis {\i)} we denote: 

\ek) = J2u,^k\i) (2) 

i 

where U is the matrix that diagonalizes H^. In addition to this coherent dynamics, the 
coupling of pigment electronic degrees of freedom to protein degrees of freedom, and the 
resultant decoherent and dissipative dynamical effects must be incorporated into the 
simulation. There is little experimental information on the vibrational dynamics of the 
protein degrees of freedom in virus-templated assemblies and hence we employ a model 
that makes minimal assumptions. The electronic- vibrational coupling is chosen to be 
linear and the vibrational degrees of freedom are modeled harmonically, resulting in the 
interaction Hamiltonian: 

where a^^i denotes the annihilation operator for mode ^ coupled to the excited state 
of pigment i, and the vibrational "bath" is described by a free harmonic Hamiltonian: 
Hy = We model the harmonic bath using an over-damped Brownian 

oscillator model [37j in the high temperature limit. The spectral density is taken to 
be Ohmic with Lorenz-Drude regularization: J{uj) = ^2 -> with reorganization energy 
A = lOOcm"^, and relaxation time 7"^ = lOOfs. The true vibrational and solvent 
dynamics of TMV-templated assemblies likely has a more complex spectral density 
than this structureless model. However, in the absence of experimental data we use this 
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simple spectral density with minimal assumptions, noting that its generic features and 
parameter values are similar to the spectral densities used for simulation of electronic 
energy transfer in natural LHCs [34j . 

In order to describe the effective dynamics of the electronic degrees of freedom 
that are of interest, we average over the harmonic bath degrees of freedom. There are 
numerous methods for doing this averaging, at various levels of approximation. In this 
work we use the modified Redfield formalism [38l [S^ to derive an effective equation 
of motion for the excitonic populations of the complex once the vibrational degrees 
of freedom have been averaged. In effect the modified Redfield approach derives a 
dynamical equation for the exciton populations: 

where R is the modified Redfield rate matrix and P{t) is a vector of excitonic 
populations of length n + 2, where n = x M is the total number of chromophores 
in the complex, and we have included two additional entries to track population lost 
via exciton loss (recombination) and excitation trapping (due to localization at the 
interface and subsequent charge separation). The exciton loss and trapping rates are 
chosen phenomenologically since there are no experimentally measured values for TMV- 
templated chromophore aggregates). See Appendix B for details. An important feature 
of this approach is that while the dynamics of off-diagonal elements of the exciton 
density matrix, which correspond to the excitonic coherences, are not explicitly followed 
in time, the components of the interaction Hamiltonian that generate them are included 
perturbatively in the calculation of the population transfer rate matrix elements, see 
Eq. (42) of [38J. 

We employ the modified Redfield formalism here for several reasons. First, it 
has been shown to be reasonably accurate across a wide range of parameter regimes, 
varying from strong to weak system-bath coupling, and is thus preferred to the standard 
Redfield approach outside of the regime of weak system-bath coupling [38j. The 
approach is particularly effective at modeling energy transfer in systems with structural 
or energetic disorder that is larger than, or comparable to, the magnitude of bath induced 
reorganization effects [39j. Second, it is an extremely efficient method of dynamical 
propagation and this efficiency will be critical for the numerically intensive optimizations 
we undertake. The primary drawback of the modified Redfield model is that it only 
allows the explicit dynamical simulation of excitonic populations and not of coherences. 
Given that the excitonic populations will nevertheless be perturbatively infiuenced by 
the coherences (as described above), this is not a significant drawback for our purposes 
since we are primarily interested in asymptotic efficiency of transport and spectral width 



of absorption, both of which can be extracted from excitonic populations (see Appendix 

0- 

Finally, we note that using the symmetries of the cylindrical stacked disk structure 
together with the dipole-dipole approximation implies that we only have to consider a 
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limited domain for the two angles that define the transition dipole of a chromophore: 
e G [0,7r/2], (j) G [-7r/2,7r/2], as weU as for the helical angle On G [-tt/M, tt/M]. The 
Hamiltonians resulting from all other choices of these angles can be replicated by a 
choice within these domains - e.g. the coupling between disks with all dipoles pointed 
up and all dipoles pointed down yields the same excitonic structure. Therefore, in order 
to improve the efficiency, we allow the angular design variables to vary within these 
restricted domains. 

3. Optimized light harvesting in cyhndrical chromophore assembUes 

3.1. Multi- objective optimization and Pareto fronts 

The multi-objective optimizations described here employ = 10 disks in the cylindrical 
aggregate, with the energy of the chromophores on any disk allowed to vary in the range 
400 — 450nm. We assume full functionalization of the TMV protein monomers and hence 
there are M = 17 chromophores per disk. All chromophores are assumed to have the 
same transition dipole strength, =3 Debye. When the two angles defining the 
orientation of the transition dipole of the chromophore, 9^(j) and the helical angle 9h^ 
are also included, this yields 13 design parameters that define the Hamiltonian of the 
electronic degrees of freedom of the chromophore assembly. We seek to optimize over this 
design parameter space in order to find the structures that are optimal for two objectives: 
transport efficiency and spectral width of absorption. The maximization of these 
objectives is performed here using a genetic optimization algorithm that evaluates the 
energy transfer efficiency and absorption spectral width for each member of its evolving 
population. [Appendix A| provides more detail on the specifics of the genetic optimization 
algorithm. Assemblies with and without disorder are considered, where the disorder can 
be both structural and energetic. For the simulation of systems with disorder, we average 
the objectives over 504 instances of disorder at each design parameter configuration. 
The disorder is introduced as independent, Gaussian distributed perturbations of the 
relevant parameters, i.e., of each chromophore's transition energy (with variance 2nm), 
transition dipole orientation angles 0, (j) (with variance 0.2 radians), and helical angle 9h 
(with variance 0.2 x 271 /M radians). Note that the introduction of disorder will break 
the equivalence of the M chromophores on a single disk. 

In the field of multi-objective optimization, a critical concept is the Pareto 
front. This is the locus or curve in objective space formed by all the solutions to 
the optimization program that optimally negotiate the trade-offs between multiple 
objectives, where this means that increasing the value of a particular objective (for 
a maximization program) , leads to the decrease in the value of one or more of the other 
objectives. The shape of the Pareto front reveals the amount of competition between 
the various objectives. In cases where there is significant competition between different 
objectives, one must choose a trade-off solution that suitably compromises between the 
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competing objectives. See Appendix A for more details on Pareto fronts and multi- 
objective optimization. Our genetic optimization is designed to converge on the Pareto 
front for the problem of light harvesting using a cylindrical antenna, with the dual 
objective functions of efficiency of excitonic transport and spectral width of absorption. 



3.2. Pareto fronts for TMV103 and TMV123 

The Pareto fronts for simultaneous optimization of spectral-width and efficiency for 
TMV103 and TMV123 are shown in Fig. [3| for both calculations with and without 
disorder, First consider the fronts for the ideal structure in the case of no disorder 
(black dots and red triangles for TMV103 and TMV123, respectively). The efficiencies 
and spectral widths achievable are very similar for both the densely packed (TMV103) 
and sparsely packed (TMV123) structures, and the shape of the trade-off curve is very 
similar too. The presence of curvature in both of these Pareto fronts indicates that there 
is a trade-off between achieving high spectral width and high efficiency. In particular, we 
see that for low spectral width absorption, both of these idealized cylindrical structures 
can achieve near unit efficiency of electronic energy transport. The similarity of these two 
fronts at all except the largest spectral widths suggests that in the absence of disorder 
there is generally little advantage achieved by dense packing and strong coupling of 
pigments; the same efficiencies and spectral widths can be achieved with dense and 
sparse packing. The exception is the regime of high spectral width, where there appears 
to be a crossover to and from a regime in which the densely packed TMV103 shows a 
better simultaneous optimization of the two objectives, consistently achieving a higher 
efficiency for a given spectral bandwidth (in the range 100-120 nm). In general, however, 
the results without disorder show that the negotiation of the tradeoff between the two 
objectives is comparable for densely and sparsely packed structures. 

With the introduction of structural and energetic disorder, several marked changes 
occur to the Pareto fronts. Firstly, both fronts collapse in the horizontal direction, 
reffecting the result that much smaller spectral widths are now achievable. There is 
also a small but noticeable amount of collapse in the vertical direction, indicating that 
somewhat smaller efficiencies are achievable with disorder. In addition, the Pareto fronts 
for TMV103 and TMV123 separate and become distinguishable. These generic changes 
can be rationalized in term of energetic and structural features, which we will do below. 
For low spectral width absorption, the densely packed TMV103 (blue squares) shows 
a slight enhancement in efficiency over the sparsely packed TMV123 (green crosses), 
although both are still very close to the optimal values obtained with the ordered 
structures. It is noteworthy that the best efficiencies achievable are comparable for 
TMV123 and TMV103, and somewhat surprisingly, not very different from the best 
efficiencies achieved in the absence of disorder. This indicates that efficiency of energy 
transfer in cylindrical structures can be robust to disorder. In contrast to this robust 
behavior of the efficiency, the achievable spectral widths are however now much smaller 
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Figure 3. Pareto fronts in the efficiency-spectral width objective space for TMV103 
and TMV123 structures with and without structural and energetic disorder. The 
statistical variation in the points with disorder is negligible (we conffimed this) at the 
level of sampling performed and therefore error bars have not been included for clarity. 
The inset is a zoom into the congested high efficiency region. 



in the presence of disorder. Here the dense packing of TMV103 is clearly beneficial 
and allows significantly greater spectral width to be attained. In general terms, the 
overall shapes of the Pareto fronts in the presence of disorder are similar for TMV103 
and TMV123, and show that a more drastic trade-off between efficiency and spectral 
bandwidth is required than in the absence of disorder. Thus, to increase spectral width 
beyond ^ 68nm 75nm) for TMV123 (TMV103) requires a significant sacrifice in 
efficiency of energy transport. 

The generic nature of the Pareto fronts obtained for the two TMV-templated 
chromophore assemblies with and without disorder indicate that the trade-off between 
efficiency and spectral width is fundamental for such cylindrical light harvesting 
antennas and independent of the presence of disorder. Beyond this result however, 
several subtle differences are apparent between the detailed behavior of the densely and 
sparsely packed systems. In particular, the densely packed TMV103 appears to be better 
able to preserve both a high efficiency and high absorption bandwidth in the presence 
of disorder. This suggests that the dense packing and ensuing quantum coherence may 
serve to provide robustness of function against sources of disorder. In the following 
section, we explore these factors in more detail by examining the parameters defining 
the structures on the Pareto fronts. 
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4. Design variables emerging from optimal structures 

To understand the Pareto fronts in Fig. [3] and to appreciate how the objectives are 
maximized and the conflicts negotiated, we now examine the performance metrics in 
parameter space. We use two complementary summaries of the structural information 
and present them in Figs. |4] and [5] for the case of no disorder, and Figs. |6] and [7| 
for the case with disorder. The flrst summary in both cases (Figs, [i] and |6]) presents 
the parameter variations for structures on the Pareto front, while the second (Figs. [5] 
and [T]) presents a detailed analysis of the excitonic and dynamical features of extremal 
structures on the Pareto fronts in Fig. [3| In the next two subsections we will interpret 
the data in these summaries, but flrst we describe in detail the content of the flgures. 

Figs. |4] and [6] summarize the values of the 13 design parameters for structures on the 
Pareto front, for TMV103 and for TMV123 without (Fig. g and with (Fig. |6|) disorder. 
The X-axis on all plots in these flgures varies over the integer parameter 5, which indexes 
structures on the Pareto fronts in Fig. [3] from left to right, i.e., starting with the most 
efl&cient (smallest spectral width) and ending with the least eflicient (largest spectral 
width). There are three panels of plots displayed in each subflgure of Figs. |4]and|6| The 
left panels, (i) show a color-coded representation of the energies of the = 10 disks for 
each structure on the corresponding Pareto front. The middle panels, (ii) show the two 
orient ational angles ^, (f) and the helical angle 6^ for each structure on the corresponding 
Pareto front. The third panels, (iii) show the ratio of coupling between neighboring 
chromophores on a single disk (intra-disk) to the coupling between closest neighboring 
chromophores on adjacent disks (inter-disk), which we deflne as A = Jintra-disk/ -^inter-disk- 
This ratio is important because it indicates the direction of strong electronic coupling 
- either within disks (large values of the ratio) or between disks and along the cylinder 
(small values of the ratio). It is important to appreciate that the electronic coupling 
(J) by itself does not dictate the rate of energy transfer between two chromophores 
(or a group of chromophores). The energy difference between them also plays a role. 
However, the value of A informs us about the direction of dominant coupling and thus 
will be used in the interpretation of results below. 

Figs. [5] and [7| present the excitonic structure and dynamics of the extremal 
structures on each of the four Pareto fronts shown in Fig. [3| The extremal structures 
are deflned as the solutions on each of the Pareto fronts that achieve the maximum 
of one of the objectives at the expense of the other: these are precisely the structures 
indexed by minimal and maximal values of s. Examining these extremal structures 
elucidates the design principles involved in maximizing the objectives. Fig. [5] examines 
extreme structures in TMV103 and TMV123 in the absence of disorder, and Fig. [7| 
examines extremes structures for TMV103 and TMV123 in the presence of disorder. 
Each subflgure in these plots is divided into two panels; the left panel shows details 
of the most eflicient structure and the right panel shows details of the structure that 
achieves the greatest spectral width. Between these two panels, in the center of the 



Design principles and fundamental trade-offs in hiomimetic light harvesting 



15 



figure we show the corresponding optimal orientations of the chromophore transition 
dipoles in the two extremal structures. Note that these orientations are shown for only 
two of the = 10 disks, the orientations on the remaining disks are identical to these. 
The structure giving the most efficient transport is on top and is to be viewed together 
with the left panel. The structure giving the greatest spectral bandwidth is below and 
is to be viewed together with the right panel. 

Considering now the left and right panels of Figs. [5]and[7| each of these shows three 
plots. In all of these plots the x-axis indexes the excitons (sorted from lowest to highest 
energy) for the corresponding extremal structure. Various properties that quantify the 
delocalization and energetics of the excitons and the resulting modified Redfield rates 
for energy transfer are shown as a function of the exciton index, as follows. 

(i) The top plot shows two measures that quantify the optical accessibility and 
delocalization of the excitons, respectively. The green line shows the magnitude 
of the excitonic dipole in units of the individual chromophore transition dipole 
magnitude (3 Debye), which constitutes a measure of the optical accessibility of the 
excitons. The blue line shows the inverse participation ratio (IPR) of the exciton 
between disks , defined as: dIPR/e = lo"^ ^2 where is the probability of exciton 
k being localized on disk n which is a sum of the probabilities of exciton k being 
localized on any of the chromophores composing disk n: p^ = X^Jt^i,^^ where 
the sum is over chromophores i that compose disk n and Ui^k are elements of the 
transformation matrix that defines the site-basis to exciton-basis transformation, 
Eq. ([2]). The disk inverse participation ratio dIPR measures the extent to which 
the excitons are delocalized over multiple disks; a dIPR value of d indicates that 
the exciton has significant amplitude over chromophores on d disks. 

(ii) The middle plot shows the energy and the spatial location of the excitons. The 
blue line denotes the energy of the exciton. The green line gives the disk number 
where most of the amplitude of the exciton is localized; i.e. arg max^ p^ for any 
exciton indexed by k. This quantity constitutes a measure of the approximate 
spatial location of each exciton. Oscillations in this localization parameter (as a 
function of k) are a heuristic 

(iii) The bottom plots show the non-zero elements of the modified Redfield rate matrix 
R describing exciton dynamics in the structure. The magnitude of each element is 
color coded. These rates are normalized by the largest rate in the matrix and so 
are all between and 1. This plot nicely illustrates the domains of connectivity 
in excitonic/energy space. Thus, if two groups of excitons are connected by large 
matrix elements in i?, then there will be significant population exchange between 
them. Individual domains are formed by groups of strongly coupled excitons that 
are only weakly coupled to other excitons. 

In the next section we will refer to these sets of plots in Figs. [5] and [7| to interpret 
and explain the details of the parametric trends shown in Figs. [4]and|6) We note that 
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(b) TMV123 without disorder 



Figure 4. Design variables for structures on the Pareto fronts without disorder^ for (a) 
TMV103 and (b) TMV123. The x-axis on ah plots shows the structure index s on the 
Pareto front. Left panels (i): color coded disk energies for the A/" = 10 disks. Central 
panels (ii): transition dipole orientation angles 9^ (j) and helical angle 9h^ with the upper 
subplot showing the variation of efficiency (green line) and spectral width (blue line) 
as a function of s. Right panels (iii): ratio of intra- to inter-disk dipole-dipole coupling, 
A. 



for the disordered cases, Figs. |6]and[7[ the parameters, structures, and excitonic details 
shown are before the introduction of disorder. That is, the physical structures, and 
consequently excitonic details, shown in these figures are perturbed by each instance of 
disorder that is sampled over. 

4.1. Optimal structures in the absence of disorder 

We first compare the optimal values of the 13 design variables for TMV103 and TMV123 
without disorder. Fig. [4j We see that the trend in variance of energies of disks is 
similar for both structures - i.e. the variance increases steadily as spectral width 
increases. Thus the incorporation of either distinct chromophores or distinct local 
energetic environments for these will be necessary to achieve large spectral width in the 
ordered structures. Furthermore the energy profiles in Figs. 4(a)| (i) and 4(b)K i) show 



that the most efficient structures have little heterogeneity in terms of pigment energies, 
whereas in the opposite extreme, the high-spectral width structures have highly varied 
pigment energies (although the bottom disk where the trap is situated always has the 
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Figure 5. Analysis of excitonic structure and dynamics for the extremal structures on 
the TMV103 and TMV123 Pareto fronts {without disorder). See main text (beginning 
of section [4]) for explanation of the plots. 



lowest energy chromophores) . The structures in the middle of the Pareto front, i.e., 
those that achieve a balance between efficiency and spectral width, have a very specific 
distribution of pigment energies that defines an energy funnel that gradually transitions 
from high energy pigments at the top to low energy pigments at the bottom of the 
cylinder. 

Turning to the choice of angular design variables, it is evident that the tangential 
angle 9 has the most infiuence on light harvesting function. For TMV103 (Fig. |4(a)| (ii)), 
the most efficient structures have 9 Q and this angle transitions to a tilted optimal 
value 0.6 as we move along the Pareto front. For TMV123, there is a similar 
transition from ^ ?^ for the most efficient structure at the left extremal of the Pareto 
front, io 9 0.4 for the extremal structure with largest spectral width, although for 
this more sparse distribution of chromophores, the corresponding transition happens 
much faster (i.e., at smaller s) and the tilted orientation is preferred for most structures 
on the Pareto front. The radial angle for all structures on the Pareto front for both 
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TMV103 and TMV123 is restricted to be in the range \(j)\ < 0.2. Therefore having 
transition dipoles that are tilted into or out of the disk plane does not seem to present 
any advantage for either efficient transport or spectral width. Within this small range 
the radial angle can fluctuate and still yield similar values of efliciency and spectral 
width. Finally, the helical angle 6^ (whose limits are smaller: — 7r/17 <9h< vr/l?) can 
also fluctuate signiflcantly without aflFecting either objective. This suggests that in the 
absence of disorder, the objectives are fairly robust to both the radial and helical angles 
(at least in the ranges \(j)\ < 0.2, — 7r/17 < Oh < 7r/17). 

To understand the behavior of Pareto structures with respect to the tangential 
angle, ^, it is useful to also look at the behavior of the ratio of intra-disk and inter-disk 
dipole couplings A, which is shown in Fig. [4](iii). It is evident that in the absence of 
disorder, this ratio drops signiflcantly as s increases and that it shows a transition from 
A > 1 to A < 1 at about the same point as 9 transitions from zero to ^ 0.6 (^ 0.4) 
for TMV103 (TMV123). This transition in A as a function of s indicates that strong 
coupling within disks is preferred for improved efliciency, while strong coupling between 
adjacent disks is preferred for enhanced spectral width. Consequently the optimal 
structures show a transition in the direction of the dominant coupling as one moves 
along the Pareto front. To conflrm that 9 is indeed responsible for this change in the 
direction of the dominant coupling within the cylindrical array, we plot the parametric 
dependence of the dipole-dipole coupling between two neighboring pigments on the 
same disk ( | Jmtra-disk | ) and two neighboring pigments on adjacent disks ( | Jinter-disk | ) , for 
TMV103 and TMV123 in Fig. |8]. We show here only the parametric dependence on the 
tangential {9) and helical {9h) angles, since these dependencies are the most relevant. As 
can be seen from this flgure, the coupling between pigments on the same disk dominates 
for most of parameter space for both TMV103 and TMV123. However, for a region 
around 9 ^ 0.6, this intra-disk coupling becomes very small and is dominated by the 
inter-disk coupling. This region where the inter-disk coupling dominates is larger for 
TMV123 because the distance between pigments on the same disk for this structure 
is larger. So we conclude that the spatial direction of electronic coupling, which is 
measured by the value of the ratio A, can be controlled by the tangential angle of the 
transition dipole orientation. 

This analysis raises the question as to what is the functional reason for transitioning 
from dominant coupling within disks to dominant coupling between disks as we progress 
along the Pareto front from most eflicient to largest spectral width, i.e., as we increase 
s? At flrst glance, this appears somewhat counter-intuitive, since one might expect that 
strong coupling between disks would lead to eflicient excitation transfer, and conversely 
that strong excitonic coupling within a disk would lead to spectral broadening. However, 
the true optimal solutions reflect a more subtle balance of the light harvesting properties 
possible in this design variable space and geometry. 

Fig. [8] shows that in both TMV103 and TMV123, due to the close proximity of 
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pigments within a disk, it is difficult for the couphng between the disks to dominate over 
the intra-disk couphng. There is only a small region around 9 ^ 0.6 where this happens 
and in this region both intra-disk and inter-disk couplings are small 10 cm~^). In 
particular, the inter-disk coupling at this angular configuration is dominated by the 
reorganization energy A = 100 cm~^ and hence the transport down the cylinder will 
be Forster-like and on slow timescales. Therefore the strategy of simply maximizing 
coupling between disks to have maximum energy transfer down the cylinder will not 
work because of the geometric dimensions in these systems. The alternative strategy 
that is converged on by the optimization of transport efficiency in our calculations, 
is to have large intra-disk coupling and smaller inter-disk coupling, while maintaining 
only a slight energy gradient. The slight energy gradient implies that neighboring disks 
contain almost resonant chromophores (recall that there is no disorder at this stage) 
and hence even the small coupling between neighboring disks creates excitons that are 
delocalized across these disks. We can confirm this strategy by examining the extremely 
efficient structures in [5} For TMV103 the delocalization of excitons by this strategy is 
corroborated by the large dIPR values and the oscillation of the disk number in which 
the excitons have primary spatial locality (both in the left panel of Fig. [5](a)). This 
network of moderately delocalized excitons with overlapping spatial locality results in a 
dense modified Redfield rate matrix (bottom panel of Fig. [5](a)) and efficient population 
transfer. The most efficient structure in TMV123 utilizes a similar strategy as TMV103, 
except that the extent of exciton delocalization across disks is a little greater due to 
the smaller energy gradient (in this case, all disks are almost the same energy). For 
small values of 5, TMV123 settles on values of 9 such that the inter-disk and intra-disk 
couplings are comparable. This creates significant energetic overlap between delocalized 
excitons that are distant and results in efficient long-range (multi-chromophoric) Forster 
transfer [40]. These observations are evidenced in the left panel of Fig. [5](b) that shows 
the significant delocalization of excitons in the most efficient TMV123 structure, and 
the dense modified Redfield matrix describing the efficient long-range transfer of energy. 

We now consider the structures with largest spectral width, shown in the right 
panels of Fig. [5j Here the strategy for optimizing spectral width in TMV103 and 
TMV123 is identical. The first element, as already noted is a varied distribution of 
transition energies: the main mechanism by which spectral width is enhanced is simply 
by having pigments that absorb at various energies. In addition to this, it is interesting 
that the angular configurations for the structures maximizing spectral width converge 
on a coupling ratio of A < 1. This results in the inter-disk coupling dominating or 
being the same order as the intra-disk coupling. To understand this, we note that 
even a small intra-disk coupling is sufficient to create delocalized excitons within disks 
(because all chromophores on the same disk are resonant), whereas a larger inter-disk 
coupling is required since chromophores on neighboring disks may be energetically 
different. Therefore, in order to create a large cluster of coupled chromophores and 
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hence excitons with large delocahzation, it is advantageous to maximize the inter-disk 
couphng while keeping the intra-disk coupling small. This is precisely the strategy 
utilized in the extremal structures with the largest spectral width, i.e., the large s 
value structures at the righthand end of the Pareto front. Indeed, the right panels of 
Fig. [5](a) and [5](b) show that this approach delocalizes excitons over multiple disks 
(as evidenced by the dIPR and disk population oscillations), which is a signature of 
electronic coupling of many chromophores over multiple disks. Such electronic coupling 
of many chromophores results in spectral broadening. It is interesting that the excitons 
that are most delocalized are the low and high energy ones. This expands the absorption 
spectrum beyond the limits set by the range of single chromophore transition energies 
- i.e. the absorption profiles extend beyond 400nm and 450nm. 



4-2. Optimal structures in the presence of disorder 

We now examine the structures and parameters of the optimal configurations on the 
Pareto fronts when energetic and structural disorder are included. The presence of 
disorder invalidates many of the optimization strategies outlined above for disorder- 
free systems. Firstly, the most striking observation for the optimal design variables of 
the disordered systems, is the presence of a strong energy funnel in both Figs. |6(a)K i) 
and |6(b)| (i). To achieve any efficiency in the presence of disorder an energy funnel 
is absolutely necessary, and in both cases (TMV103 and TMV123), a shallow energy 
funnel amounting to a gradual decrease in pigment energies as the disks approach the 
bottom layer is more effective at generating efficient transport. With the exception of 
these most efficient structures, the energetic profiles for the remaining structures on 
the Pareto front for both TMV103 and TMV123 are remarkably similar. The optimal 
design strategy for the non-extremal structures appears to be to have discrete decreases 
in energy every couple of disks, leading to a more step-like energy gradient than the 
smooth gradient preferred by the most efficient structures. 

Turning to the choice of angular design parameters, consider the case of TMV123 
first (Fig. |6(b)| (ii)). Unlike the disorder-free case, where a variety of angle choices 
(especially for (j) and 9^) yielded similar objectives, in the presence of disorder a 
few optimal configurations emerge. In particular, the most efficient structures have 
9 = 9^ ^ 0. At the other extreme, the structures with the largest spectral width have 
9 ^ 7r/2 and 9h ^ 0. The radial angle is constrained to be 10 1 < 0.25 and is close to 
zero for nearly all structures on the Pareto front. As with the disorder- free case, we find 
that the tangential angle has the biggest infiuence on the objectives. The radial and 
helical angles have only a slight infiuence, as long as they are constrained to be ^ 0. For 
disordered TMV103 (Fig. 6(a)| (ii)), the tangential angle again has the biggest infiuence 



and it transitions from 9 ^ 0.4 to ^ 7r/2 as 5 increases. The other angles fiuctuate 
along the Pareto front of TMV103 more than for TMV123, but we cannot discern any 
pattern in their variation along the Pareto front. Finally, Figs. |6(a)| (iii) and 6(b)| (iii) 
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(a) TMV103 with disorder 




(b) TMV123 with disorder 



Figure 6. Design variables for structures on the Pareto fronts with disorder^ for (a) 
TMV103 and (b) TMV123. The x-axis on ah plots shows the structure index s on the 
Pareto front. Left panels (i): color coded disk energies for the = 10 disks. Central 
panels (ii): transition dipole orientation angles 6>, (j) and helical angle Oh^ with the upper 
subplot showing the variation of efficiency (green line) and spectral width (blue line) 
as a function of s. Right panels (iii): ratio of intra- to inter-disk dipole-dipole coupling, 



show that there is a clear preference for having dominant coupling along the cylinder 
(low A) to optimize efficiency, while dominant coupling within disks is preferred in order 
to optimize spectral width. In fact the most efficient structures maximize the inter-disk 
coupling. For TMV103 this is achieved around 6 ^ 0.4, 6>/, = and for TMV123 it is 
achieved around ^ 0, ^/^ (see Fig. |8]). And similarly, the structures with maximum 
spectral with maximize the intra-disk coupling, which occurs at 7r/2 as seen from 
Fig. 11 

We can gain more insight into the parameter landscape by examining the extreme 
structures on the disorder Pareto fronts. As the analysis below will show, the inability to 
construct large delocalized excitonic states on any timescale in the presence of significant 
static disorder dominates the design choices made in optimizing light harvesting. The 
left panels of Figs. |7(a)| and |7(b) show that the most efficient structures for TMV103 
and TMV123 have similar excitonic structure, showing several distinct features. Firstly, 
the excitons are mostly localized on a single disk. There is a slight derealization, 
especially in TMV103, as evidenced by the dIPR and oscillations in disk population 
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Figure 7. Analysis of excitonic structure and dynamics for the extremal structures 
on the TMV103 and TMV123 Pareto fronts {with disorder). See main text (beginning 
of section [4]) for explanation of the plots. 



but it is minimal compared to the disorder- free structures. This is a result of the 
strict energy gradient, with each disk having pigments with a different transition energy 
and consequently little delocalization across disks. The addition of static disorder, and 
also dynamic disorder induced by protein fluctuations, will of course further localize 
the excitons. The structure of the modifled Redfleld rate matrix for the most eflicient 
TMV103 and TMV123 structures (bottom plots of left panels Fig. [?[ is seen to be 
fairly sparse with small rates between densely coupled domains. The domains are 
localized on disks and there is larger transfer rates between excitons within a domain 
than between domains. This localization of excitons into domains and inter-domain 
dominated transport is reminiscent of multi-chromophoric Forster transfer [40j , and we 
can expect this to provide accurate description of the energy transfer in these structures 
optimized for efliciency in the presence of disorder. The strategy here will be to have 
slightly delocalized excitons on disks that are as strongly coupled as possible to excitons 
on neighboring disks due to maximized inter-disk electronic coupling (at 6 ^ 0.4, 9h = 
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for TMV103 and at ^ 0, ^/^ for TMV123), and a shallow energy gradient which 
ensures that neighboring disks contain chromophores that are nearly resonant despite 
the disorder. 




e (rad) e (rad) 



(a) TMV103 (b) TMV123 

Figure 8. Dipole-dipole coupling between two neighboring pigments on one disk 
(brown-tan surface) and two neighboring pigments on adjacent disks (red- yellow 
surface) as a function of the tangential transition dipole orientation {0) and the inter- 
disk helical {Oh) angles (the radial transition dipole orientation (j) is fixed here at zero), 
(a) TMV103, (b) TMV123. The primary difference between these two structres is the 
density of pigments within a disk. The distance between pigments on adjacent disks for 
both TMV103 and TMV123 is 20l. The distance between pigments on the same disk 
is - (27r/17) X 25 = 9.24 Afor TMV103 and - (27r/17) x 40 = 14.78 Afor TMV123. 

Now we turn to the extremal structures that maximize spectral width in the 
presence of disorder. These are illustrated with the lower structures shown in the central 
panels of Fig. [7| There is little obvious difference between these two extremal structures 
for TMV103 and for TMV123. Both consist of disks that are J-aggregate-like, with 
dipoles aligned head-to-tail. The energy of each disk is different, except for the first 
two and last two which have almost identical energies. This creates a series of almost 
independent J-aggregates (i.e. bright excitons are the lowest energy ones) that absorb 
at a range of energies [§J The reason for this form of the optimal structure for achieving 
maximal spectral width can be understood by examining Fig. [8| which reveals that 
the maximum possible excitonic coupling between any two chromophores (for TMV103 
and TMV123) is an intra-disk coupling that results from 6 = 7r/2, i.e., from dipoles 
aligned head-to-tail within a disk. This strong coupling is advantageous for competing 

§ Note that the excitonic details shown in Fig. [7|are for the ideal design parameters. Each instance of 
disorder will perturb these and result in a perturbed version of these excitonic details. In particular, the 
exciton delocalization for the maximum spectral width structures (for TMV103 and TMV123) seems 
large from the dIPR values shown in the right panels of Fig. [T] However, any instance of disorder will 
break the symmetries of the ideal design and result in more localized excitons. 
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against disorder and maintaining excitonic coupling in its presence. It is not, however, 
advantageous for efficient transport down the cyhnder since the strong couphng is intra- 
disk rather than inter-disk, but it is nevertheless advantageous for spectral broadening. 
We can rationalize this finding with a simple model for the effect of disorder on coupling 
of any two chromophores. Thus, given excited states of two chromophores with energies 
£"1, E2 and electronic coupling J, the resulting excitonic states have energies 

Er + E2 ^ - + 4 J2 

e± ± — . 

^2 2 

Assuming that neither are dark states [i.e. have zero net dipole strength), the absorption 
peaks of the coupled system will be centered around e±. In presence of finite coupling 
J, these energies can be different from Ei and £"2, causing spectral broadening due 
to excitonic coupling. In the presence of disorder, — £2! can be large, in which 
case e± will only be significantly different from the original transition energies {Ei 
and E2) ii J is comparably large. This is the reason for the J-aggregation seen in the 
maximum spectral width structures; this mechanism achieves the strongest coupling and 
thus enables broadening even in the presence of disorder. The spectral widths achievable 
with disorder for TMV103 are slightly larger than those achievable for TMV123 because 
the former system can attain greater intra-disk coupling strengths due to the greater 
proximity of pigments within a given disk (see Fig. [s]). In summary, in the presence 
of disorder, maximal spectral width is achieved by having each disk effectively act as 
its own antenna with strong intra-disk coupling present to provide spectral broadening 
and with almost every disk composed of pigments with distinct transition energies. 
Interestingly, this optimized structure bears resemblance to multi-junction or tandem 
solar cells [4T1|42], which are also designed to increase the range of wavelengths of utilized 
photons. 

It is interesting to note that the strategies found here for optimizing spectral 
width with and without disorder are only effective in their respective cases. That is, 
the disorder strategy of creating J-aggregate-like maximally coupled disks would not 
be effective in the absence of disorder. This is because of the well known property 
of J-aggregates that they possess red-shifted, narrow lines of absorption because the 
structural symmetry renders most states optically inaccessible {i.e. results in dark 
states). Therefore in the absence of disorder the absorption profile would consist of 
thermally broadened peaks around the J-aggregate absorption lines, which is unlikely 
to produce a very wide spectrum. However in the presence of disorder, the complete J- 
aggregate symmetry is broken and the majority of states are no longer dark states, 
allowing wide bandwidth absorption and effective spectral broadening at the same 
time (note that the maximum disorder should be comparable in magnitude to the 
maximum achievable electronic coupling for this strategy to be effective). Disorder 
is thus essential to the success of this strategy. Similarly, the disorder-free strategy that 
is optimal for maximizing spectral width is to use the spectral broadening resulting from 
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excitonic coupling across many chromophores. This strategy is in turn ineffective in the 
presence of disorder, since this would suppress such couplings and resulting broadening, 
as seen from the above equation for a two chromophore system. These observations 
underscore the important point that the optimal strategies for light harvesting can 
change dramatically depending on the amount of disorder present in the system. 

5. Design principles for engineering light harvesting antennas 

Given the understanding we have gained of the structures on the Pareto front and their 
relative performance, we can now extract several design principles for optimal light 
harvesting using cylindrical chromophore assemblies. 

(i) A trade-off between efficiency and spectral width exists for cylindrical light 
harvesting antennas. There are limits to maximizing both of these objectives, 
especially in the presence of disorder. 

(ii) Optimal structures for efficiency or spectral width change significantly depending 
on whether energetic and structural disorder is present or not. In the absence of 
disorder, many choices of the design variables yield similar values of the objective 
functions. However, when disorder is present only a few choices of the design 
variables are truly optimal - i.e., the design variables for neighboring structures on 
the Pareto front are fairly similar. 

(iii) An energy funnel is essential regardless of the density of pigments for such structures 
with a single dedicated sink for excitons. We see the emergence of an energy funnel 
for both TMV103 and TMV123. A shallow energy funnel is the most effective for 
transport. 

(iv) The attachment angle that dictates the tangential angle of the dominant dipole {9) 
is the most critical orientation degree of control, followed by the helical angle 9h. 
The radial angle should be kept as close to zero as possible since all structures 
on the Pareto front retain this property. 

(v) An increased density of pigments does not aid in increasing efficiency or spectral 
width in the absence of disorder. However, in the presence of disorder, increased 
density is advantageous since it can lead to stronger electronic couplings which can 
be used to combat the deleterious effects of disorder. 

(vi) Simply having access to strong electronic coupling is not sufficient to overcome the 
deleterious effects of disorder. The direction of this coupling has to be controlled and 
used effectively. In the case of TMV103 and TMV123 the strongest couplings are 
intra-disk simply because of the dimensions. These couplings only help to maintain 
spectral width in the presence of disorder and do not directly generate efficient 
energy transfer (although other strategies help to generate efficient transfer). 



Design principles and fundamental trade-offs in hiomimetic light harvesting 26 

We note that in the present system with a fuU complement of M = 17 chromophores 
per disk, there does not appear to be any apparent benefit from having a helical twist to 
the chroomophore arrangement along the cylinder. This may simply reflect the relatively 
close packing of chromophores in this system and the ability of the other 12 parameters 
to simultaneously optimize the energy efliciency and spectral width. However it may also 
reflect the lack of any chiral constraint element at the chromophore level. More extensive 
calculations with smaller values of M and more reflned models of chromophores are 
needed in order to determine whether additional structural constraints such as a center 
of chirality on the chromophores are necessary for achieving helical optimal structures 
for these design objectives. 

We can now ask how many of these design principles are manifest in natural 
photosynthetic light harvesting antennae. Such observations are unavoidably speculative 
since we do not understand the structure-function relationships in natural antennae 
sufliciently well at this stage to make general statements. Nevertheless, it is still 
interesting to compare these emergent design principles to features and motifs seen 
in natural systems. Firstly, consider chlorosomes, the main light harvesting antennas of 
green sulfur bacteria, which are cylindrical antennas consisting of very densely packed 
BChl d, or e molecules [SI [4]. Given the above results, it is plausible that the dense 
packing found in chlorosomes has the primary function of increasing the range of usable 
photons, an important feature for green sulfur bacteria, which typically live in severely 
energy limited environments. It is diflicult to speculate on the impact of structure on 
energy transfer efliciency because it is unknown at present how the photo-excitation 
couples out of the chlorosome complex and hence we cannot model energy extraction 
simply as trapping sites at the base of the cylinder as done above. One could also 
speculate that the small energy gradient seen in PS-II [43] is consistent with the above 
observation above that small energy gradients are beneflcial for maximum efliciency of 
transport in disordered systems. 

6. Conclusion 

Inspired by the cylindrical assemblies of chromophores that can be synthesized by TMV- 
templated assembly, we have examined the landscape of light harvesting performance 
achievable by such structures. In particular, we have studied the trade-oflFs involved 
in optimizing over multiple objectives relevant to light harvesting. In addition to 
identifying a fundamental trade-off between optimizing energy transfer efliciency and 
bandwidth of absorption, our calculations have allowed the development of several 
design principles to guide the design and construction of such cylindrical assemblies 
of chromophores. We flnd that the presence of disorder drastically eflFects the optimal 
structural forms and therefore it is essential to characterize the amount of structural 
and energetic disorder present in typical TMV-templated chromophore assemblies. 
Experiments to-date on protein-templated systems are insuflicient to characterize this 
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degree of disorder, and hence this is a critical task for future experiments. The 
design principles established here provide guidelines for a program of quantum-informed 
rational design (QuIRD) for biomimetic light harvesting systems. 

Our study also lends insight into the structure of biological photosynthetic light 
harvesting complexes. The multi-objective optimization study clarifies the potential role 
of strong electronic couplings enabled by high densities of chromophores. In realistic 
systems with disorder, our study suggests that strong electronic couplings are somewhat 
beneficial for maintaining efficiency of energy transfer, but more importantly, they 
are essential for increasing the spectral width of absorption profiles. Therefore the 
strong electronic couplings recently observed in several biological LHCs - e.g. Refs. 
[441 I45l [461 WJ\ ^ may not only aid excitation transport in such systems, but also 
simultaneously benefit optical measures of light harvesting performance such as spectral 
width of absorption. 

This first multi-objective optimization of structural design for biomimetic light 
harvesting systems suggests several interesting directions for further work. The first 
avenue for future work is the consideration of measures of performance other than 
absorption bandwidth. Although efficiency of energy transfer is likely to be desirable 
in all applications, other measures of performance in the optical domain could be 
considered. For example, sensitivity to photons of a particular wavelength could be 
relevant for the design of sensor technologies. As another avenue of future work, we 
plan to carry out a multi-objective optimization study of the helical structural that 
is templated using TMV [33]. The helical assembly could be more stable than the 
stacked disk assembly considered in this work, however features such as an energy 
gradient could be more difficult to implement in the helical structures. Finally, we plan 
to also perform multi-optimization studies of TMV-templated chromophore assemblies 
that are restricted to varying the experimental degrees of freedom that are currently 
tunable. That is, the design variables will be restricted to the ones that are currently 
experimentally accessible, so that the optimized designs could be constructed and 
verified in the immediate future. 
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Appendix A. Multi-objective optimization 

The optimization of multiple objectives is a commonly encountered problem in 
engineering and sciences. Given a set of control variables, x = (xi, X2, Xat), the 
task of multi-objective optimization, in this example a maximization, is: 

max [/i(f),/2(f), .../m(^)] (A.l) 

X 

where fi{x) are M objectives (cost functions). One could, and usually does, also 
have equality or inequality constraints on the variables. For convenience we notate 
/ = /2(^), .../m(^)]- This optimization over multiple objectives has multiple non- 

degenerate solutions unlike a single objective optimization which has a single solution 
(or multiple degenerate ones). The set of solutions are referred to as Pareto points 
and are the set of points in objective space such that improving any one objective 
can only be done at the expense of another. Specifically, / is a Pareto point, or is 
Pareto optimal, if there does not exist another feasible objective vector such that 
fl > fi Vi G {1, 2, M}, and /j > fj for at least one j G {1, 2, M}. 

The set of Pareto points in objective space is also called the Pareto front. For each 
Pareto point in objective space, /*, there is a corresponding Pareto point in control 
variable space, x* and is specified by the control variables that achieve /*. 

One approach for solving multi-objective optimizations is to combine the multiple 
objectives into one objective, e.g. F{x) = ^^fLi^ifii^) with a G M is a linear 
combination of objectives. Then a conventional single objective optimization solution 
method is used to solve this problem. This strategy is acceptable when the acceptable 
combination of objectives is known a priori. However, in many problems one cannot 
know ahead of time how to combine the multiple objectives into one. In addition, 
knowledge of the Pareto front for the full multi-objective problem is valuable for 
understanding the types of trade-offs involved in the optimization problem. In this 
work, where we are particularly interested in characterizing the trade-offs involved in 
light harvesting, the full multi-objective optimization ins critical. 

There are a variety of methods for solving multi-objective optimizations [15j. In 
this work we employ an evolutionary algorithm which has the advantage that it requires 
few assumptions and a priori knowledge of the optimization landscape. We utilize the 
non-dominated sorting genetic algorithm II (NGSA-II) by Deb et al [49] implemented 
in the C++ optimization toolbox written by Sastry |50j. The specific genetic algorithm 
parameters we use are: 

(i) Population size: 100 

(ii) Number of generations evolved: 100 

(iii) Proportion of population replaced in each generation: 0.7 

(iv) Crossover probability (with simulated binary crossover): 0.85 

(v) Mutation rate (with polynomial mutation): 0.3 
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These parameter values were converged upon by experimentation and determining 
what combination of population size, mutation rate and crossover rate allowed good 
exploration of the optimization landscape. For each optimization, multiple runs were 
executed (there were several restarts), each from a different randomly chosen starting 
population and all the resulting populations were collected. This distribution of starting 
points added another degree of randomization to ensure that the large optimization 
landscape was sampled reasonably well. 



Appendix B. Modified Redfield theory 

The equations of modified Redfield theory specify rates of exciton population transfer, 
and define a rate equation for the exciton population dynamics: 

'-^ - RPit) (B.l) 

where R is the modified Redfield rate matrix |38], and P{t) is a vector of exciton 
populations. In this work we only consider exciton dynamics in the single exciton 
manifold since this is the most relevant at low to moderate solar irradiance. As derived 
in Ref. [38j, the rate for population transfer from exciton k to k' is: 

poo 

i?,,,, = 2Re / rfrF* (r)A,(r)A,,,,(r) (B.2) 

with 



Ffc/(r) = exp(-i(£;^, - Afc/)r - g^k' ^k'k'{^)) 
Ak{T) = exip{-i{E^ + Xk)r - gkkM^)) 

Nk,k'{^) = {gk'k,kk'{^) — [gk^k,kk{^) — gk'k,k'k'{^) — ^^i^'k.k'k'] 

X [gkk'Mk{r) - gkk'M'k'{r) - 2i\kk'^k'k']) e2(^fcfc,fe/feKr)+^A,,,,,,r) 

Here, = E^ — Xk is the 0—0 exciton transition energy ^ i.e. E^ is the exciton transition 
energy and Xk is the exciton reorganization energy, defined by A^^^^ = ^^^^ \Un,k\^Xn with 
Un^k the exciton-site basis transfer coefficients and the reorganization energy of site 
n. Similarly, Ac,^^^^^ = '}2n=i ^n,a'^n,/3U^^^Un,s^m and the exciton lineshape function is 
defined as ga(3,js{t) = Yln=i Un,aUn,/3U^^^Un^sgn{t) where gn{t) is the single chromophore 
(site) phonon-induced lineshape function. In this work we assume all chromophores 
have the same lineshape function, and it derives from a high temperature over-damped 
Brownian oscillator model of the phonons with spectral density J{uj) = ^2 • In this 
case. 




(B.3) 
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The integrals defining the modified Redfield rates are time-consuming to perform 
numericaUy for an arbitrary hneshape function, but in the case of the model above, the 
integral defining the modified Redfield rates becomes of the form: 

poo 

Rk^k' oc 2Re / dr e^^^-^^^~^^ (B.4) 

where q are fc, k' dependent, but time-independent, complex coefl&cients. Integrals of 
this form can be performed analytically using a continued fraction expansion, as shown 
by Takagahara et al in the Appendix of Ref. [48j. We use this continued fraction 
solution to calculate the modified Redfield rates efficiently. 

Now we detail how to extract a measure of spectral width of absorption and 
efficiency of transport from the modified Redfield model of exciton dynamics. 

Appendix B.0.1. Spectral width An expression for the absorption spectrum in terms of 
the modified Redfield exciton transfer rates is [51j: 

A{uj) =ujJ2 d^Rej r dt exp (i{uj - ujk)t - I UnM'dit) - ^) } (B.5) 

k n=l 

where the k sum is over all excitons. dk and ujk = Ej^/h are the magnitude of the 
transition dipole and the transition frequency of exciton fc, respectively. g{t)^ as detailed 
above, is the phonon induced lineshape function, taken to be that of an over-damped 
Brownian oscillator and the same for all chromophores. Finally, is the inverse lifetime 
of exciton fc, given by a sum of outgoing modified Redfield rates: jk = — ^k'^k ^k'k'kk- 
We construct the absorption spectrum according to this expression and normalize 
it to have maximum 1. Then the width of the absorption spectrum is defined as the 
sum of frequency intervals for which A{uj) > 0.1. 



Appendix B.0.2. Efficiency Consider the modified Redfield rate equation, Eq. (B.l), 
which describes the dynamics of the exciton populations. In this work, we use 

P{t) = (]9i(t),]92(t), ...,px(t),Ptrap(^)'^^loss(^))^' where the first K elements are the 
populations of the K excitons at time t (in the main text we consider cylindrical TMV 
assemblies with disks and M chromophores per disk, and therefore K = N x M). 
The remaining two elements in this population vector are the population in the trap 
(i.e. the fraction of excitons that have undergone charge separation), Ptrap(^)' 
the population lost to recombination, Pioss(^)- elements of R that represent rates 
of transfer between excitonic populations are prescribed by modified Redfield theory as 
above, and the rates of transfer from exciton k to the trap (^/^^trap) ^^^^ (^fc^loss) 
channels are defined as: 



|2 



^fc^trap ~ 7trap WnM 
nedisk N 

^fc^loss =7lossd^ (B.6) 
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The motivation behind that k trap rate is that exciton k has probabihty 
^nedisk N being locahzed at one of the chromophores on the bottom most 

disk (disk A^) where the trapping (charge separation) is assumed to take place. Hence 
the total trapping rate is this probability multiplied by the bare trapping rate, 7trap- 
this work we choose l/7trap ^ ^ps. The motivation behind the k loss rates is that the 
rate of radiative recombination of an exciton is equal to a bare rate of recombination 
7jQgg multiplied by the dipole strength of the exciton. This ensures that bright and 
dark excitons recombine at different rates. In this work we choose the bare rate of 
recombination as VtIoss ^ 

We define the efficiency of transport as the asymptotic fraction of population in 
the trap, rj = Ptrap(^ ~^ Using a Laplace transform solution of the rate equation 
Eq. (B.l) results in the following expression for this asymptotic population: 

P{t ^ oc) = lim(5l - R)-^P{0) (B.7) 

where P{0) is the initial population distribution. We approximate this asymptotic 
solution by taking a small value of s in the above expression. We have confirmed that 
our choice of s does not infiuence the solution, especially the value oirj = Ptrap(^ ~^ 

Finally, our choice of P(0), the initial exciton distribution is dictated by the size 
of each exciton's dipole. That is, Pk{0) = d^/A/", where A/" is a normalized so that 
J2k=iPk(^) ^ ^- T^^^^ initial state specifies that the probability that an exciton is 
initially populated is proportional to its transition dipole strength. The initial trap and 
loss populations are set to zero. 



